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Z . Then the error in this computed solution is 
E — z — z 


ABSTRACT 


We study the error in the derivatives of an unknown function. We construct 
the discretized problem. The local truncation and global errors are discussed. 
The solution of discretized problem is constructed. The analytical and 
discretized solutions are compared. The two solution graphs are described by 
using MATLAB software. 

KEYWORDS: Differential Equations, MATLABHeat Equation 

1. INTRODUCTION 

1.1. Measuring Errors 

In order to discuss the accuracy of a numerical solution, it is necessary to 
choose a manner of measuring that error. It may seem obvious what is meant by 
the error, but as we will see there are often many different ways to measure the 
error which can sometimes gives quite different impressions as to the accuracy 
of an approximate solution. 

1.1.1. Errors in a Scalar Value 

First we consider a problem in which the answer is a single value z^R. 
Consider, for example, the scalar ODE 
u’(t)=f(u(t)),u((>)=n (i-U 

and suppose we are trying to compute the solution at some particular time T, so 
z=u[ T). Denote the computed solution by 


( 1 . 2 ) 


1.1.2. Absolute Error 

A natural measure of this error would be the absolute value 
of E, 

|i?| = |z-z| (1.3) 

This is called the absolute error in the approximation. 

1.1.3. Relative Error 

The error defined by 


t -» 0 ■ This is slightly stronger than the previous statement, 
and means that f(j) decays to zero faster than g(r)- If 
f (f) zz o(g(r)) then /(r) = 0(g(r)) though the converse may 
not be true. Saying that f(j) = o(l) simply means that the 
/(r)^Oasr^O. 


Examples 1.1.1 

2r 3 = 0(r 2 )as T - 


2 3 1 

0 , since = 2r<\ for all r < -. 

2 2 


is called relative error. 


1.1.4. “Big-oh” and “little-oh” notation 

In discussing the rate of convergence of a numerical method 
we use the notation Of: p ), the so-called "big-oh" notation. If 
f (t) and z( T ) are two functions of T then we say that 

/rU£r»«~ 0. 


2r 3 = o(t ) as r — » 0, since 2r —» 0 for all T < — ■ 

2 

r 3 t 5 

sin(r) = 0(t) as T —» 0 ,since s in(r) = r - — + — + ... < z for all 

T > 0 . 

sin(r) = r + o(t) as T —» 0, since ( sm ( r )~ r ) = 2 ) . 


1.1.5. Taylor Expansion 

Each of the function values of u can be expanded in a Taylor 
series about the point x, as e.g., 


If there is some constant C such that 


m 


g(j) 


<C 


for all T 


sufficiently small, or equivalently, if we can bound 
\f{r)\<C\g(r)\ for all T sufficiently small. 


u(x + t) = u(x) + m'(x) + — r 2 u\x) + — t 3 u'"(x) +0(t 4 ) 
2 6 


u(x -t) = u(x) - m'(x) + —t 2 u"(x) - — r 3 u m (x) + 0(r 4 ) 
2 6 


(1.4) 

(1.5) 


It is also sometimes convenient to use the "little-oh" 


notation/( r ) = 0(g(r)) as r -» 0. This means that 


m 


gw 

/ 


^as 


1.2. Finite Difference Approximations 

Our goal is to approximate solutions to differential equation, 
i.e. to find a function (or some discrete approximation to this 
function) which satisfies a given relationship between 
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various of its derivatives on some given region of space and/ 
or time along with some boundary conditions along the edges 
of this domain. A finite difference method proceeds by 
replacing the derivatives in the differential equations by finite 
difference approximations. This gives a large algebraic 
system of equations to be solve in place of the differential 
equation, something that is easily solved on a computer. 

We first consider the more basic question of how we can 
approximate the derivatives of a known function by finite 
difference formulas based only on values of the function itself 
at discrete points. Besides providing a basis for the later 
development of finite difference methods for solving 
differential equations, this allows us to investigate several 
key concepts such as the order of accuracy of an 
approximation in the simplest possible setting. 


= u\x ) + t 2 u\x ) + 0(t a ) 

= u" (x) + 0(t 2 ) 

1.2.3. Higher Order Derivatives 

Finite difference approximations to higher order derivatives 
can be obtained. 

S + S 2 u(x) = \[u(x + 2t)-3u(x + t) + 3u(x)-u(x-t)] (1-10) 

T 

= u m (x) + ^m"\x) + 0(z 2 ) 

= u"'(x) + 0(t ) 


Let w(x] represent a function of one variable that will always 
be assumed to be smooth, defined bounded function over an 
interval containing a particular point of interest x, 


The first one equation (1.10] is un-centered and first order 
accurate: 

S 0 S + S_u(x) = —^ [u(x + It) - 2 u(x + t) + 2 u(x - t) - u(x - 2 r)\ (!■ H) 


1.2.1. First Derivatives 

Suppose we want to approximate u\x) by a finite difference 

approximation based only on values of u at a finite number of 
points near x. One choice would be to use 

e . u(x + t)-u(x) 

d + u(x) = - 

T 

= u'(x) + ^m"(x) + 0(T 2 ) t 1 - 6 ) 

for some small values of T. It is known as forward 
difference approximation. 

Another one-sided approximation would be 

e x u(x)-u(x-t) 
o_u(x) = —^ 1 

T 


= u’\x) + -t 1 u""\x) + 0(t 4 ) 

4 

= u\x) + 0(t 2 ) 

This second equation (1.12] is second order accurate. 

2. Comparison Of Analytical and Discretized Solution 
Of Heat Equation 

2.1 Solutions for the Heat Equation 

2.1.1. Finite Difference Method 

We will derive a finite difference approximation of the 
following initial boundary value problem: 

U t = U xx for x e (0,l),f > 0 

u(0,t) = u(l,t) = 0 for ? > 0 ( 2 . 1 ) 
u(x, 0) = f{x) for xe (0,1) 


= u(x)-^nt\x) + 0(r 2 ) 


(1.7) 


for two different points.lt is known as backward difference 
approximation. 

Another possibility is to use the centered difference 
approximation 

u(x+t)-u(x- 
2 T 



1 

2 


[S + u(x) + S_u(x)] 


Let n > 0 be a given integer, and define the grid spacing in the 
x-direction by 


The grid points in the x-direction are given by 
Xj = j r for j=0,l,...,n+l. Similarly, we define t m =mh for 

integers m>0> where h = At denotes the time step. Then, we 
let 

V J denote an approximation of u(x 7 >t m ]. We have the 

following 

approximations 


= u'(x) + -t 2 u'" + 0(t 4 ) 
6 


U-8) U ^)- U(X ’ t + h) ~ U(X ’ l) +0(h) 

h 


( 2 . 2 ) 


1.2.2. Second Order Derivatives 

The standard second order centered approximation is 
given by 

/ x u(x-t)-2u(x) + u(x + t) m 

o z/(x) =--- 


and 

, x u(x-T,t)-2u(x,t) + u(x + T,t) f'J 

= —--- — + 0(t 2 ) (2.3) 

T 

These approximations motivate the following scheme: 
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v;. 1 -2v;+v; +1 


for j=l,...,n, m >o 


C2.4) 


By using the boundary conditions of (2.1), we have 
Vq = 0 and v" +1 =0 , for all m > o _ 


The scheme is initialized by 


V? =/(*,) 


J J , for j=l,...,n. 

h 

r = — 

Let T . Then the scheme can be rewritten in a more 
convenient form 


v™ +1 = rv™ + (1 - 2r)v™ + rv‘ 


™>0 


(2.5) 


When the scheme is written in this form, we observe that the 
values on the time level tm+1 are computed using only the 
values on the previous time level tm and we have to solve a 
tridiagonal system of linear equations. 


2.1.2 Approximate Solution 

The first step in our discretized problem is to derive a family 
of particular solutions of the following problem: 


vr-v; y;- 1 - 2 v;+v; + , 

^ for j=l,...,n, m ^ 0 (2.6) 

with the boundary conditions 
v m = 0 v m = 0 

0 and n+l , forall m -°. (2.7) 

The initial data will be taken into account later. We seek 
particular solutions of the form 

v m — YT 

j j m for j=l,...,n, m ^ 0 (2.8) 

Here X is a vector of n components, independent of m, while 

{T } 

1 is a sequence of real numbers. By inserting (2.8) 

into (2.6), we get 

X j T m + i -*jT m _ Xj_ x T m -2XT m +X J+1 T m 
h r 2 


T m +1 - (1 h ^) T m f or m > 0. (2.12) 

Some iterations of (2.12) 


T m+l =(\-h M )T m 
= (\-h J u) 2 T m _ l 


= (l-hjuT +l T 0 
= (1 -hju) m+l 


Clearly indicate that the solution is 

T m=( l ~ h M) m f or m > 0. (2.13) 

This fact is easily verified by induction on m. Next we turn 
our attention to the problem (2.9) with boundary condition 
(2.11). In fact this is equivalent to the eigenvalue problem. 

Hence, we obtain that then eigenvalues ^ are gi ven 

by 


2 2 

= — —-cos(tc^t) 

T T 

2 

= — (1 - cos(ktzt)) 




for 


k = 1, 


■>n 


(2.14) 


and the corresponding eigenvectors 

X k =(X Ktli X Ka ,...,X K JeR", K=l _ n 


have components given by 
X K,j = sin( Kjg r) 


Since we are looking only for nonzero solutions, we assume 
XT ^0 

that 7 m , and thus we obtain 


T -T 

1 m +1 _ m 

hT 


X 


7-1 


-2X j + X j+l 


t 2 X , 


X Kj =sm(K7g Xj ) ^ . n 

It can be easily verified that 

AX k j = 4 x ^+A x .j ~\ x *,m 

T T T 


The left-hand side only depends on m and the right-hand side 
only depends on j. Consequently, both expressions must be 

equal to a common constant, say^ - ^, and we get the 
following two difference equations: 


X, 


- 2 X. +X . 


= ~VXj 


, for j=l,...,n, 


(2.9) 


~ T rn 

h 


= -MT m 


for m - 9. 


( 2 . 10 ) 


: ~sm(far(j - l)r) + ^sin (wgr) - ^sin(*7zr(y + \)r) 

T T T 


1 2 1 
= —r-sin(A:^r - K7tt)) + — si n{K7rjr) -r-sin(A:^r + kttt) 


sin (K7qx) cos {kuz) - cos (K7ijr) sin(/c^r) - 2 sin (jC7gr) 
+ sin (jC7gr) cos(/c;zt) + cos (jC7gr) sin(/c;zT) 


We also derive from the boundary condition (2.7) that 

X o _ X n+ 1 - 0 (2.11) 

We first consider the equation (2.10). We define T0=1 and 
consider the difference equation 


= — sin(/c^’r)(l - cos(atzt)) 
t 

= \ sin(*oy r)(2 sin 2 (^)) 
r 2 
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( 2 . 20 ) 


4 . 2 KKT 

- — sin (-)sm(A7Z7r) 


v m . 

Hence, we obtain particular solutions K,J of the form 


X K = ( KJC ) 2 for K = 1,2, 


= (1 -hju K ) m sin (k7dc.) 


(2.15) 


with the eigenfunctions 
X k (x) = SW(K7BC) for K = 1,2,. 


T&— T 

Further, the solution of * is given by 

*• = 1,2 . 


( 2 . 21 ) 

T K (t) = e^ for 

( 2 . 22 ) 


We have derived a family of particular solutions 1 with 

V m . 

values K ’ J at the grid point (xj,tm). Next, we observe that 
any linear combination of particular solutions 

n 

K=l , (^ is scalar] 


is also a solution of (2.6) and (2.7). Finally, we determine the 
coefficients { } by using the initial condition 

for )=!,...* 

Since v *" — ^ at t=0, we want to determine {} such that 

Y,r K x K ,j= f(*j) 


, for j=l,...,n. 


(2.16) 


Hence, it follows from 

n 

f(- x i ) = X sin (ktdCj) 

K=\ 


X sjj( x j)) h =YjYM nK7Dc j’ s ' ms7Dc i 


that 

Yk =2 T^f(Xj)X K 


M 


' for 


2.1.3. Exact Solution 

To find a solution of the initial-boundary value problem (2.1), 
assume that 

u(x,t) = X(x)T(t) (2.17) 


Finally, we use the superposition principle to a solution to the 
initial-boundary value problem (2.1), and we get 

00 

u(x,t) = Y i c K T K (t)X K (x) 

K =1 


u(x,t) = ^c K e (K7r) 1 sin (k7zx) 

K =1 


The unknown coefficients can be determined from the initial 
condition such that 

QQ 

/W = 2X sin^) 

*=> (2.23) 


The solution of the continuous problem is given by 

u(x,t ) = y^c K e~ Aj sin(ATzx) 

(2.24) 


where ^ ( K7r ) fourier coefficient 

i 

c K = 2 J / (x) sin(/r;Dc)<ix 

0 (2.25) 


2.2. Comparison of Analytical and Discretized Solution 
2.2.1. We want to compare this analytical solution with 
the discretized solution given by 


Vj = sin (KfiXj) 

K =1 


(2.26) 


where 



(2.27) 


Using boundary conditions we get 
X(0)=X(1) = 0 


If we insert the (2.17) in the equation (2.1), we have 
X\x) _ ftt) 


x(x) m 


(2.18) 


Now we observe that the left hand side is a function of x , 
while the right hand side just depends on t. Hence, both have 
to be equal to the same constant This yields the 

eigenvalue problem for x 

X" = -AX 1 X(0) = X(l) = 0 (2.19) 


and 


n 

Y K = 2r V f(Xj ) sin (kot, ) 

M f()r K-\,...,n ( 2 . 2 8) 

In order to compare the analytical and discretized solution at 

u m = u (x. t ) 

a grid point (xj,tm), we define j v J ’ , i.e, 

oo 

U 1 = ^j C K e ~ Ktm sin(A : 7 aj) 


(2.29) 


Our aim is to prove 


V, 


u. 


2 

Nontrivial solutions only exist for special values of *■ 

They are so-called the eigenfunctions. In this special case we 
have the eigenvalues 


under appropriate conditions on the mesh parameters T and 
h. To avoid technicalities, we consider a fixed grid point 

(xj,tm) where ~ ^ for t >0 independent of the mesh 
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parameters. Furthermore, we assume that the initial function Since the sine part here is identical, it remains to compare the 

f is smooth and satisfies the boundary conditions, i.e. f(0) Fourier coefficients c * and and the time-dependent terms 

=f(l) =0. Finally we assume that the mesh parameters h and 


u m , v m 


T are sufficiently small. 

u m 

In order to compare J and J , we note that 

oo 

«(*,■ ^ m ) = Z c * e_V " sin(A :m j ) 

K=\ 

= sin(m.)+ sin(/cflx y .) 

x*=1 /f=n+l 

Here we want to show that 

oo 

sin (/ctdc.)^ 0 


(2.30) 


Since f is smooth, it is also bounded and then the Fourier i n fact we have 

c i 

coefficients K are bounded for all K . 

Obviously, we have 

|sin(Am y .)| < 1 


and 

oo 

Z c v e ~ V " ■&m(Knx j ) 


OO 

- SKI| e ' v '| sin (^)| 

K=n+\ . . 

—y— 

<i 


< max cj > e 


-(ktt) t m 


< c 


c 


ik'k. 


t,„ >t 


7 V+l / 2, V- 

') +r'J 


+ ... 


= cl 

:Cl 




l + e Kt +e lKt +... 


~ 0, for large value of n. 

Since we have verified (2.30) it follow that 

n 

u™ ~ AJm sin(^cc y -) 

-i (2.31) 

Now we want to compare the finite sums (2.26) and (2.31): 

n 

vj = Z Y* C 1 - y sin (^y ) 

K =1 

Motivated by the derivation of the solutions, we try to 
compare the two sums term wise. Thus we keep K fixed, and 
we want to compare 

c K e~ Kt " sa^taacj) and y K (l-hju K ) m sm(»a,)_ 


e K " and ^ . 

2.2.2. Comparison of Fourier coefficient c * and 
coefficient 

We start by considering the Fourier coefficients, and note that 
Y* is a good approximation of c * because 

n 

2 ^f(Xj)sm(K7!Xj) 

2-1 

is the trapezoidal-rule approximation of 

i 

2 J / (x) sin(A :7Dc)dx 
o 

2 J / (x) sin (K7ix)dx - 2 / (x y ) sin(ftm y ) 




M 


2 + / (x y . )]sin(/rTzx 7 ) - 2r^ / (x y .) sin(KT 2 x 7 .) 

z j =i 7=1 


r Z F _ r ) “ /( x 2 )]sm(Km,) 


7=1 


=-z 

7=1 L 


/(*,) - /(*, ) 


sin(Km.) 


Z 

7=1 L 




|sin(*m y .)| 


^Z 


-T 2 f(x j ) + — f n (x j ) + ... 


) , for f sufficiently smooth. 


-Xj m 


2.2.3 Comparison of the Terms and e 

We will compare the term ^ ~~ approximates the term 
, we simplify the problem a bit by choosing a fixed time 


e Kn 


h = — 


tm, say tm=l and we assume that 2 . As a consequence, 
we want to compare the terms 


<* K =e K 

and 

P K =(1-^/0 


(2.32) 


(2.33) 


Since both and ^ are very small for large values of K , it 
is sufficient to compare them for small K . In order to 

compare a « and P* for small values of ^*, we start by 
recalling that 

sin(y) = y + 0(y 3 ) 


Thus, we get 


2sin | I = 2 


h 

ten A— - 


(^Tl ?! 2 (^) 5 f ^ 2 


3! 


5! 
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= 2(»t) 2 | 


3! 5! 


h t 

IJ =U,(x J,t m ) +-u„(Xj,t u xx (x J ,t m )-—u xxxx (x j ,t m )-... 

h t 2 

= uA X j>tJ-U xx (Xj,tJ + -U t t(Xj,tJ- — U x ^(Xj,t m )-... 


« 2{k7t) 2 — 

2 , for h sufficiently small 
= K 2 7T 2 h 

By using these facts, we derive 





( 

r 

l*T 

1 - 2 sin 2 

Kn i 

V 

v v 



A 


h t 

= 2 “# ( x 7 + - 

= 0(r\h) 

The finite difference scheme (2.6) is consistent of order 

( 2 , 1 ). 

Example 2.3.1 

Solve the 1VP 

u t =u xx f or x e(0,1) t> o 

u(()j) = u(\.l) = 0 f Q| . / > () (2.34) 

u(x,0)=f(x) for xe(0,l) 


P K =(1 -K 2 7t 2 h) 


= 1 + -(-icVA) + 
h 


{-K 2 7T 2 h) 2 + ... + (-K 2 7T 2 h) h 


= \-K 2 7T 2 + ] -^-(-K 2 7T 2 ) 2 +... + (-K 2 7T 2 h) h 

~ o~ Klnl 

~ e , for h sufficiently small 

= e~ K 

P K ~ a K 

-X K t m 

This shows that also the time dependent term e is well 
approximated by its discrete ^ ~ [ 

2.3 Consistency 
Lemma 2.3.1. 

The finite difference scheme (2.6) is consistent of order (2,1). 
Proof: 

Local truncation error of the finite difference scheme (2.6) is 
given by 

ur-< i/ x 


The exact solution is given by 
a = -e~ n 1 sin 7ix 


Approximate solution and exact solution are 



Solutions 


7,7 m m+1 „, m ( m ) 

hi. = Uj -Uj + u j+1 ) 


h 2 h 3 

= u(Xj,t m ) + hu t (Xj,t m ) + —u tt (Xj,t m ) + —u m (Xj,t m ) +... 


-u{Xj,t m )- — 


u(X:,t )~ TU (x -,t m ) H- U (X/,t ) - U ( X j,t ) 

j 2 J 6 J 

r 4 

+- U xxxx (Xj,t m ) -+...-2u(Xj,t m ) +u(Xj,t m ) + TU x (Xj,t m ) 

t 2 t 3 r 4 

- U (x .,t ) H- U (x .,t ) H- U (x .,t ) + ... 

2 J 6 xxx 7 m 24 xxxx y m 


MATLAB codes for Figure 2.1 
N=100 ; x=zeros (N,l); 
g= zeros (N,l); 
v=zeros (N,N); 

%delta x=h, delta t=p 

%f (x)=cos (2*pi*x); 
H=1/(N+1); 

P=1/(N+1); 


= hu,(Xj,t m ) +—u tt (Xj,t m ) + ... 


z 2 u (x J )h- u (x J )- 

xx V i ’ m / xxxxv /’ 


for j = 1:N 
x(j)=j*h; 

end 

%for m=l:N 
%t(m)=m*p; 
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terml=g(k]. *v(k,:); 
suml=suml+terml; 
s =suml; 
end 


%end 
%gamma k 
for k=l:N 

sum=0 ; 
for j=l:N 

term=2*h*sin ( 2*pi*x (j ))*sin (k*pi*x (j )); 

sum=sum+term; 

end 

g ( k ]=sum ; 

end 

%vk 

for k=l : N 
for j = 1 : N 

v (k, j )=p*(4/h A 2)*(sin(k*pi*h/2]] A 2)%in(k*pi*x(j)); 
end 

v(k,:); 
end 

sum 1 =0; 
for k= 1 : N 


%%%Approximate solution 

plot (x, s ,' - ') ; hold on 

%%%%%%Exact solution 

Plot (x,-exp (-*pi A 2*0.01 )*sin(pi*x)); 

3. CONCLUSION 

The aim of this research paper describe the Errors between 
Analytical solution and Discretized solution of Differential 
equations. 

4. REFERENCES 

[1] R. J. LeVeque, Finite difference methods for differential 
equations, Washington, 2006. 

[2] J. W. Thomas, Numerical partial differential equation: 
Finite difference methods, Springer-Verlag, New York, 
1995. 


@ IJTSRD | Unique Paper ID - IJTSRD27937 | Volume - 3 | Issue - 5 | July - August 2019 


Page 2272 




